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Ph ■ Abstract 

For a family of second-order elliptic operators with rapidly oscillating periodic coef- 
ficients, we study the asymptotic behavior of the Green and Neumann functions, using 
Dirichlet and Neumann correctors. As a result we obtain asymptotic expansions of 
Oh I Poisson kernels and the Dirichlet-to-Neumann maps as well as near optimal conver- 

■ gence rates in W l ' p for solutions with Dirichlet or Neumann boundary conditions. 

1 Introduction 



The main purpose of this paper is to study the asymptotic behavior of the Green and 
Neumann functions for a family of elliptic operators with rapidly oscillating coefficients. 
More precisely, consider 



O 
^1" 



C £ = -div(A(x/e) V) = ° 



° y \e) dx 



j J 



6>0 (1.1) 



dxi 

O . 

(the summation convention is used throughout the paper). We will assume that A(y) 



{djf{y)) with 1 < i,j < d and 1 < a, (3 < m is real and satisfies the ellipticity condition 

m 2 < agivmf < -ICI 2 for y G R d and £ = (C) e M dm , (1.2) 

/i 

where \i > 0, and the periodicity condition 

A(y + z) = A(y) for y G R d and z G Z d . (1.3) 

We will also impose the smoothness condition 

\A(x) - A(y)\ < t\x - y\ x for some A G (0, 1] and r > 0. (1.4) 

Let G e (x,y) = (G^(x,y)) and N £ (x,y) = [N^(x,y)) denote the Green and Neumann 
functions respectively, for C e in a bounded domain Q, with pole at y. We are interested in 
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the asymptotic behavior, as e — > 0, of G E (x,y), N e (x,y), V x G e (x,y) and V x N e (x,y) as well 
as V x V y G e (x,y) an V x V y N £ (x,y). We shall use G (x,y) and N (x,y) to denote the Green 
and Neumann functions respectively, for the homogenized (effective) operator £ i n ^- 

Let Pj = Xj(0, . . . , 1, . . . ) with 1 in the (5 th position for 1 < j < d and 1 < (3 < m. 
To state our main results, we need to introduce the matrix of Dirichlet correctors $^ • = 
. . . , in n, defined by 

^(Sj,.) = in n and = Pf on dfi, (1.5) 

as well as the matrix of Neumann correctors j = (^j, • • • , ) in ^, defined by 

£ £ (^.) = in ft and _(^.) = onffl. (1.6) 

Here djdv e denotes the conormal derivative associated with C e for e > 0. 
The following are the main results of the paper. 

Theorem 1.1. Let C £ = — div (A(x / e)^) with the matrix A(y) satisfying conditions U.ty) . 
U.3\) and Ji.^| ). Then for any x, y G Q, 

Cf 

\G e (x,y)-G (x,y)\ < (1.7) 

\x — y\ 

if Q is a bounded C 1,1 domain, and 

if Q is a bounded C 2,v domain for some rj G (0, 1), where C depends only on d, m, fi, r, A 
and Q. 



Theorem 1.2. Suppose that A(y) satisfies the same conditions as in Theorem li.il Also 

hf(y) = a fi 



assume that A* = A, i.e., a"f(y) = a?°(y) for 1 < i,j < d and 1 < a, < m. Then for any 



x,y G tt, 

\N E (x,y)-N {x,y)\ < 'V?' ' ~ J ( L9 ) 



Celn 


\e- X \ 


X 


-y\ 


+ 2] 




x — 


y\ 


d-i 



if Q is a bounded C 1 ' 1 domain, where C depends only on d, m, fi, r, A and Q. Moreover, if 
Q is a bounded C 2 ' v domain for some r\ G (0, 1), 

\^{Nf( x ,v)} - £{«J3<*)} ' £m*«)}\ < C '^%llt 2] (i-xo) 

for any x,y G £1 and t G (0, 1), where M = diam(ft) and Ct depends only on d, m, fi, t, X, 
t and Q. 

A few remarks are in order. 
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Remark 1.3. In the case of a scalar equation (m = 1), the estimate (11. 7p holds for bounded 
measurable coefficients satisfying (II. 2p and (11.31) (see Theorem 13.31) . 

Remark 1.4. The matrix of Dirichlet correctors ■) was introduced in [2] to establish 
the boundary Lipschitz estimates for solutions with Dirichlet conditions, while the matrix 
of Neumann correctors (^f was introduced in [21] to establish the same estimates for 
solutions with Neumann boundary conditions. It is known that ||$^ ■ — P^\\l°°(q.) — Ce and 
IIVC, II L°°(n) + || V\I / f J -||L°°(n) < C Under the condition ^^(xo) = Pj(xo) for some x$ G Q, 
we also have ||^ - Pf\\ L <*>(n) < Ce\w[e~ l M + 2] (see Propositions GOO and Q. 

Remark 1.5. Estimates ( jl.7j) and ( 11.9j) in Theorems 11.11 and 11.21 allow us to establish 0(e) 
estimates for \\u e — Uo\\lp(q) (1 < P < oo) for solutions with Dirichlet or Neumann boundary 
conditions (see Theorems 13.41 and I4.5p . More importantly, estimates ( 11. 8p and (ll.lOp yield 
near optimal convergence rates in W l,p (Q) for any 1 < p < oo. In fact, let C £ (u £ ) = F in fl 
and u £ = on dQ. Then 

IK -no- «• - Pf}^\\ w ?V) < C pz{ He~ l M + 2}}^ l \\F\\ LP{n] . (1.11) 
In the case of Neumann boundary conditions we obtain 



u 



e -u - - pP}-^JL\\ wl , pm < Ct^WFWvM (1.12) 



for any t G (0,1), where C £ (u £ ) = F in Q, |^ = on dQ and J Q F = J m u £ = 0. See 

subsections 3.2 and 4.2 for details. Let w £ = u £ — u — exj(x/e)^-, where {Xj{y)) denotes 
the matrix of correctors for Ci in M d . Estimates ( II. lip and ( I1.12p should be compared to 
the well known 0(e 1 ^ 2 ) estimate: ||w £ ||#i(^) = 0(e 1 / 2 ) (see e.g. [6|), and to the following 
estimate, 

1/2 



+ |y \Vw £ (x)\ 2 dist(x,dn)dx^ < Ce\\F\\ L 2 {n)j (1.13) 

proved in [201 Theorems 3.4 and 5.2]. Due to the presence of a boundary layer, both the 
0{e 1/2 ) estimate and ([USD are more or less sharp. The Dirichlet and Neumann correctors 
are introduced precisely to deal with boundary layer phenomena in periodic homogenization. 

Remark 1.6. Our approach to Theorems 11.11 and 11.21 also leads to asymptotic estimates 
of V x V y G £ (x, y) and V x V y N £ (x, y) (see subsection 3.3 and Remark l4.9p . As a result we 
obtain asymptotic expansions for {d / dxi){C £ )^ 1 {d / dxj) and A £ , the Dirichlet-to-Neumann 
map associated with C £ . 

The asymptotic expansion of the fundamental solutions as well as the heat kernels for 
C £ in R d has been studied, using the method of Bloch waves; see e.g. [2S1 IUJ and their 
references (also see [23] for results obtained by the method of G- convergence). In the presence 
of boundary, the Bloch representation is no longer available. In a series of papers [21 HI 
[3], M. Avellaneda and F. Lin introduced the compactness methods, which originated in 
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the regularity theory in the calculus of variations and minimal surfaces, to homogenization 
problems. In particular, they established an asymptotic expansion for V y G £ (x, y) in [I], using 
Dirichlet correctors. As a result, it was proved in |3] that if Q is C l,r) for some rj G (0, 1), 

P £ (x, y) = P (x, y)oJ £ (y) + Re(x, y) for x E Vt and y G d£l, (1.14) 

where P £ (x,y) (e > 0) denotes the Poisson kernel for C £ in Q, 

d 

"e(y) = g^{®tk(y)} • Mv)h(y) ■ ^(i/KO/Ki^A) (1.15) 

and the remainder R e (x,y) satisfies 

limsup [\R £ (x, y)\ : x G E and y G dVt] = (1-16) 

for any compact subset E of Q (the results are stated for the case m — 1; however, the 
argument in |4j works equally well for elliptic systems). In f 1 1 . 1 5 j) we have used $* k to 
denote the Dirichlet correctors for £*, the adjoint of C £ . Also, h{y) = (aijni(y)nj(y))^ 1 
and {chij) is the (constant) coefficient matrix of C Q . The expansion (I1.14p was used in |1] to 
identify the limit, as e — > 0, of solutions to a problem of exact boundary controllability for 
the wave operator ^ + C e . 

Our Theorem 11.11 gives a much more refined estimate of R £ (x,y) in (jl,14p (under the 
stronger condition dQ G C 2,v ). Indeed, it follows from the estimate (II. 8p that 



Celnfe- 1 ! 


x - 


- y\ 


+ 2] 




x - 


-y\ 


d 



R £ (x,y)\< — 1 — -p for any x G fi, y G dVt. (1.17) 



Besides its applications to boundary control problems, estimate (11.171) may also be used to 
investigate the Dirichlet problem 

C £ {u £ ) = in Q and u e (x) — f(x, x/e) on dfl, (1-18) 

where f(x,y) is 1-periodic in y. The Dirichlet problem f 1 1 . 1 8 [) arises natually in the study 
of boundary layer phenomena and higher-order convergence in periodic homogenization (see 
e.g. [251 [U H51 EH] and their references). Let v £ be the solution to 

£q{ v e) — in Q and v e = f(x, x/e)oj e (x) on dfl, (1-19) 

where u £ is given by (I1.15p . It follows from the estimate (11.17)) that 

\K - v £ \\ LP{n) = 0((e[\n(l/e)] 2 ) 1/p ) for any 1 < p < oo 

(see Theorem I3.9p . This effectively reduces the asymptotic problem (ll.ISp to the study of 
convergence properties of u £ on dQ, under various geometric conditions on Q. This line of 
research will be developed in a future work. 

We now describe the main ideas in the proof of Theorems 11.11 and 11.21 The basic tools 
in our approach are representation formulas by Green and Neumann functions, uniform 
estimates for Green functions established in [2], 

\V x G £ (x,y)\ + \V y G £ (x,y)\<C\x-y\ 1 - d and | V x V y G £ (x, y)\ < C\x - y\- d , (1.20) 
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and the same estimates obtained in |20j for Neumann functions N e (x, y). Let D r = D(xq, r) = 
B(x ,r) fl Q and A r = A(x ,r) = B(x ,r) D dfl for some x G Q and < r < r . First, to 
establish ( II. 7} , we will show that if p > d, 

IK - Uo|U°°(Di) < C / |u e - u Q \ dx + Ce\\Vu Q \\ L oo {Di) + C p e\\V 2 u \\ LP ( D4 ) } (1.21) 
Jd 4 

where C £ {u £ ) = C (u ) in _D 4 and u £ = u on A 4 . This is done by considering w £ (x) = 
u e (x) —uo(x) —ex(x/e)Vuo and using the Green representation formula and the observation 

that C £ (w £ ) = e-j^- {b%jk (x/e) gf 2 g° fc J , where bijk(y) is a bounded periodic function. Estimate 

( II. 7p follows from ( ll.2ip by a more or less standard argument (see subsection 3.1). Next, we 
show in subsection 3.2 that 

C f 

\\Vu £ - (V$ £ )(Vm )|| L o O(Di . ) < — - / \u £ - uq\ dx + Car' 1 \\ Vu \\ L °°(D 4r ) 

r J D Ar (1.22J 

+ Celn[e-V + 2]|| VV)IU~(D 4r ) + Cer"|| V 2 u \\ c o,v {Dir ), 

if C £ (u £ ) = Co(uq) in D^ r and u £ = uq on A4 r . Estimate ( II. 8p follows easily from fl 1 . 2 2 [) by 
taking u E (x) = G £ (x,yo) and u (x) = Go(x,yo)- By repeating the argument, estimate (ll.22p 
also gives an asymptotic expansion for V x V y G £ (x, y) (see Theorem 13. lip . To prove fl 1 . 2 2 j) . 
we let 

3 ' dxj 

and represent w £ in D r , using the Green function in D, where D is a C 2,ri domain such that 
D 3r C D C D 4r . 

Although a bit more complicated, the proof of Theorem 11.21 follows the same line of 
argument as Theorem ll.il In subsections 4.1 and 4.2 we establish boundary L°° and Lipschitz 
estimates similar to f 1 1 . 2 1 j) and (ll.22p for u £ and uq satisfying C £ {u £ ) = £o(uo) in D^ r and 
fir = on ^ 4? " ^ ne resu ^ s rer y on the uniform LP and Neumann function estimates 
obtained in (231 IZI] under the additional symmetry condition A* = A. 

The rest of the paper is organized as follows. Section 2 contains some basic formulas 
and estimates which are more or less known. The case of Dirichlet boundary conditions is 
treated in Section 3, while Section 4 is devoted to the case of Neumann boundary conditions. 
In Section 5 we prove two inequalities, which are used in subsection 4.3 and are of interest 
in their own right, for the Dirichlet-to- Neumann map A . 



w £ = u £ {x) - u {x) - - if }^ (1.23) 



2 Preliminaries 

Let C £ = -div(A(ar/e)V) with A(y) = (a°f(y)) satisfying Q-(0]). Let X (v) = (xfiv)) 
denote the matrix of correctors for C\ in R d , where Xj(y) = (Xj(y)i • • ■ ■> xT(v)) is defined 
by the following cell problem: 



A(x?) = -MP? 



in 



(2.1) 
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for each 1 < j < d and 1 < < m. Here Y = [0, l) d ~ M. d /Z d and P?(y) = yj(0, . . . , 1, . . . , 0) 

with 1 in the (3 th position. The homogenized operator is given by Cq = — div(v4V), where 
and 

dy. (2.2) 



A = (a^ p ) and 



a ij 



Y 



It is known that the constant matrix A is positive definite with an ellipticity constant de- 
pending only on d, m and /i (see [6]). 
Let 

btf(y) = a;? - off (y) - ^(y)^- k (xf)- (2.3) 

Since J Y bfdy = and ^-(Kf) = by Q and (j2TT]) . there exists F$ E H\Y) such that 

b «=wS F ^ and ^ = - F ^ (2 - 4) 

Remark 2.1. To see flUD, one solves Af°f = bf in Y with e ff^K) and / y = 0, 

and let 



k ^ y> dy k dy t 



(see e.g. [20]). Note that if A(y) is Holder continuous, then Vx and hence 6°f are Holder 
continuous. It follows that VF is Holder continuous. In particular, ||x||ci(y) + ll-^llc* 1 ^) i s 
bounded by a constant depending only on d, m, fi, A and r. In the case of the scalar equation 
(m = 1) with bounded measurable coefficients, the corrector \ is Holder continuous by the De 
Giorgi -Nash estimates. This, together with Cacciopoli's inequality and Holder's inequality, 
implies that there exist t > and C > 0, depending only on d and /x, such that 

/ |Vx| dy < Cr d - 1+t for x e Y and < r < 1. 

J B(x,r) 

In view of (12.31) we obtain 



/ K(y)\dy < Cr d - 1+t for x 6 Y and < r < 1. 



(2.5) 



J(«,r) 

Since A/y = % in K and j Y fydy = 0, 

IIV^IIl-^ < CWVfijWw+Csup [ T^%[dy < C{d,fj), (2.6) 

x& Jy F — 2/| 

where we have used (12.51) to estimate the last integral in (12.61) . It follows that ||-Ffcjj||oo < 
C(d,fi). 

The following proposition plays an important role in this paper. We mention that formula 
(12. 8p with Vfj(x) = P^(x)+eXj{x/e) is known and may be used to show that \\u e — ^o||l 2 (C7) < 
C^H Wo 1 1 #2^), where C £ {u £ ) = C (uo) in Q and u £ = u on dfl (see e.g. |20j). The proof of 
our main results on the first-order derivatives of Green and Neumann functions will rely on 
(12. 81) with the matrices of Dirichlet and Neumann correctors respectively in the place of the 
functions Vf 4 . 



Proposition 2.2. Suppose that u £ G H 1 (Q),uq G H 2 (Vt) and C £ (u £ ) = £o( M o) ^ n ^- 



: (x)=« e (x)-ti (x)-{^.(x)- J Rf(x)} 



(2.7) 



w/jere = (F^, . . . , V™f) G H^Q) and C e (Vf d ) = in Q for each 1 < j < d and 
1 < j3 < m. Then 



(W =e^{[i?Z(x/ e )] 



ax 



+ 



dxjdxk 

d 



d 2 ul \ 
dxjdxk J 



(2.8) 



<9xj L ' £ 

where 5^ = 1 if (3 = 'j, and zero otherwise. 
Proof. Note that 



dW 



. dxidxk 



4 {x/e) 



dxj 



a/3 



'.r\ du 1 aj3 



x\ du 







'x\ d 



n V : ) dxj ~ y V : ) dxj V : ) dxj 1 



9uq 
dx k 



a/3 



d 2 U? 



dxkdxj 



This, together with C £ (u £ ) = C (uq), gives 



^a/3 a/3 / / \ 

d 



1 5uJ 
dxj 



+ {C £ (V £ ^ k -P^)} 



a <9Mq 

dx k 



dxj 



V, 



e.k 



- x k 8^)- 



2„,7 



d u 



d 



dxidxk 
d 2 ul 



Since 
we obtain 



dxkdxj 

£s{v;[ k - P2) = -C £ {P2) = C £ {e X l(x/e)}, 



Vt k (x)-x k 5^-^(x/e) 



dxj 



8 2 ul 
dxidx k 



(2.9) 



+ 



_d_ 

dxi 
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d 2 ^ 



dx k dxj J ' 



where bff(y) is defined by (12.31) . In view of (12. 4p . we may re- write the first term in the right 



hand side of (12. 9p as 



d 8 



dxj I dxi 



kij 



J I fcy dx k dxj 



The formula ( 12. 8 p now follows. □ 
The next proposition will be used to handle the Neumann boundary condition (cf. |21j). 
Proposition 2.3. Let w £ be given by \2. 7| j. Suppose that ^{Vf k } = ^{-Pf }■ Then 

^/(^{y^-p^}-— — . (2.io) 



dv £ ) \dv £ J \du ) v e ' dx k dxj 

Proof. Note that 

<9w e \ a (du £ \ a f du \ a al3 d r a p-, dul 



d 2 ui 



(2.11) 



Since -^{Vf k } = ^{P* }, the third term in the right hand side of (l2TTT|) equals - 



+ 



fe) a This gives (l2TTUj) . □ 

The following two propositions provide the properties of the Dirichlet and Neumann 
correctors needed in this paper. 

Proposition 2.4. Let C £ = — div(A(x/e) V) with A(y) satisfying ( CQ| ) ; ( 11.51) and fll.^p . Lei 

(^ej) denote the matrix of Dirichlet correctors for C £ in a C 1,r) domain Q. Then 



V$J(x)| < C, |$^-(x) - P/(x)| < C£ (2.12) 



and 



|V{$f»-/f(x)-exJ(a;/e)}| <Cmin|l,^| (2.13) 
for any i6fl, where S(x) = dist(a:, dQ) and C depends only d, m, fi, X, t and Q. 

Proof. The first estimate in ( 12. 12ft follows from the Lipschitz estimate in [2]. To see the 
second estimate, let u £ {x) = Q^Ax) —P?{x) —eXj{x/e). Then C £ {u £ ) = in £7 and u £ {x) = 
—eXj{x/e) for x G dVl. It again follows from [2] that ||w e ||.z,°°(n) < C||ii e ||£°°(an) < Ce. Hence, 
11$^- - P/||£cx>(fi) < Ce + ||oo < Ce- Finally, note that ||Vu £ ||L°°(n) < G. Also, by the 
interior estimate in [2] and ||u e ||i«>(n) < C £ > one obtains |Vn £ (x)| < Ce[5(x)] _1 . This gives 
the estimate (l2~T3l) . □ 
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Proposition 2.5. Let C £ = — div(A(x/e)V) A(w) satisfying M.S\) , M.3\) , and t/ie 

symmetry condition A* = A. Let (^ £ j) denote the matrix of Neumann correctors for C 



in 



a C 1 '^ domain Q. Suppose ^ (xq) = if (xo) for some Xq G O. JTien 



|W?(z)| < C, |*fj(ar) - if (x)| < Celn^M + 2] (2.14) 



and 



|V{*?» - if (x) - exf (x/e) } | < Cmin |l, ^ J (2.15) 

for any x G Q, where M = diam(fi) and C depends only d, m, fi, X, r and Q. 

Proof. The estimate (12. 15ft as well as the first estimate in (12. 14ft was proved in |2Tj . To 
prove the second estimate in (12. 1411 . we let H 1 ^- = \lf ^(x) — if (a;) — eXj{x/e). Since 
| Viijf -(x) | < C min (l, £r[o~(x)] _1 ) , by the Fundamental Theorem of Calculus, we may deduce 
that \H% d (x) - H^(y)\ < Celn^M + 2] for any x,y G SI. Since \H^(x )\ = e\x?(x )\, we 
obtain \H^(x)\ < Celnfe^M + 2] for any x G f2. This gives the desired estimate. □ 

3 Asymptotic behavior of Green functions 

The goal of this section is to prove Theorem 11.11 We also establish several convergence 
theorems for solutions with Dirichlet boundary conditions. 

Let C £ = -div(-A(ac/e)V) with A(y) satisfying conditions ( Ojl and flO). Let 

G £ (x,y) denote the matrix of Green's functions for C £ in a bounded domain Q. It follows 
from [2] that if f2 is C 1 ^ for some 77 G (0, 1), 



C 

\G £ (x,y)\<- -^min., 

\ x — y\ I \ x ~ y\ \ x ~ y\ \ x ~ y\ 



6(x) 6(y) 6(x)S{y) \ 

' ; i . .i j i _. . .i j i _. io f 



iV.G.Cx.y)! < j C , d r min(l,^- 

f — y| I — y| 

V„G e x,y < | pmin l,^, 

| a; — w| d 1 { \x — y\ 

\V x V y G £ (x,y)\< ( 



(3.1) 



\x-y 



d 



for any x, y G 0, where £(x) = dist(x, 90). These estimates, which are well known for second- 
order elliptic operators with constant coefficients, play an essential role in our approach to 
Theorem 11.11 



3.1 L°° estimates 

In this subsection we give the proof of the estimate (11.71) . As a corollary of (ll.7p . we also 
establish an 0(e) estimate for \\u e — wo||lp(q) for any p > 1 (see Theorem I3.4[) . Throughout 
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the subsection we will assume that A(y) satisfies conditions (11.2p - (11.3p and in the case m > 1, 
A(y) is Holder continuous. Let 



D r = D r (x , r) = B(x , r) fl Q and A r = A(x , r) = B(x , r) fl dVL 
for some Xq G fi and < r < ro- 

Lemma 3.1. Assume that Q is Lipschitz if m = 1, and C 1,r) for some r\ G (0, 1) if m > 1. 
T/ien 

|l-(a.) < C||/||L o(A 3r ) + — / \u £ \dx, (3.2) 



where C £ (u £ ) = in D% r and u £ = f on A 



3r ■ 



Proof. By rescaling we may assume that r = 1. The estimate is well known in the case 
m = 1 and follows from the maximum principle and De Giorgi -Nash estimate. If m > 1 
and / = 0, estimate (13.21) follows directly from [21 Lemma 12]. To treat the general case, 
consider C £ (v £ ) = in D with the Dirichlet data v e — f on dD n dQ and v £ = on <9-D \ <9f2, 
where D is a C 1,ri domain such that D 2 G D G D 3 . Note that D^e ||z,°°(o 2 ) < C||/||l°°(A3) by 
[21 Theorem 3], and u £ — v £ may be handled by [21 Lemma 12], as before. □ 

The next lemma provides a boundary L°° estimate. 

Lemma 3.2. Assume that Q satisfies the same assumption as in Lemma \3.1\ Let u £ G 
H 1 ^^) and uq G W 2 ' p (D iT ) for some d < p < oo. Suppose that 

C £ {u £ ) = Cq(uq) in _D 4r and u £ = u on A^ r . 

Then, 

C f 

IK - uo|U°°(zv) < — / |« e -wo| cfe + Vwo|U=°(D 4r ) 

r ^D 4 r (3.3) 
+ CpEr^P ||V 2 Mo||LP(D 4r )- 

Proof. Note that if C £ {u £ ) = F, then C £ / r (v) = Fi, where v(x) = r~ 2 u £ (rx) and Fi(x) = 
F(rx). Thus, by rescaling, it suffices to consider the case r = 1. To this end we choose 
a domain D, which is C 1,v for m > 1 and Lipschitz for m — 1, such that D 3 G D G D4. 
Consider ^ 

w e = u £ - m - (x/e) = + wf 1 in D, 

where 

£ 6 (twW) = C £ (w e ) inD and G Hq(D) (3.4) 

and 

£ e ( w ( 2 )) = in 5 and = w £ on 9D. (3.5) 
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Since = w £ = —ex{x / e)VuQ on A3 and ||x||oo — G, it follows from Lemma [3.11 that 



ki 2) IU«(Di) < C£||Vno|| L oc (A3) +C 

Jd 3 



W\dx 



< CVIIVwoIIl^As) + C / \w £ \dx + C \w { p\dx 

Jd 3 Jd 3 



<C j \u £ — u \ dx + C*£:|| Vmo||l° c (D 3 ) + C||^e IU° 
./.D3 

This gives 

IK - Uo\\l^(D!) < C j \u E - u \ dx + C£\\Vu \\l°°{d 3 ) + C\\wP\\l°°(d 3 )- (3-6) 
Jd 3 

To estimate w^P in D$, we use the Green function representation 

wi 1 \x)= ld e (x,y)C e (w e )(y)dy, 
Jd 

where G e (x,y) denotes the matrix of Green functions for C £ in D. Using (j2.8p with Vfj = 
Pj(x) + EXj(x/e), we obtain 

)2 



= -e / — {<5 £ (x, y)} • + a y (y/e) Xfc (y/e)] • 5 dy, 

in <jy% oyjOyk 



where we have suppressed the superscripts for notational simplicity. Note that by Remark 
12.11 ||i*}jfc||oo — C. It follows that 



D 

i/ P > 



\wP(x)\ <Ce i\V y G £ (x,y)\\V 2 u (y)\dy 

(3.7) 



<Ce\\V 2 u a \\ LP{Di) ^JjV y G £ {x iy )\ p 'dy^ 

< C p 8 1| V 2 M || Lp(D 4 ) 



if p > d, where we have used Holder's inequality and || V y G £ (x, — Cp- This, together 

with (13.61) . completes the proof. □ 

We are now ready to prove the first estimate in Theorem 11.11 Note that in the scalar 
case m — 1, no smoothness condition on A(y) is needed in the following theorem. 

Theorem 3.3. Suppose that A(y) satisfies conditions M.fy) - in~!3\) . Ifm>l, we also assume 
that A(y) is Holder continuous. Let Q be a bounded C 1 ' 1 domain. Then 

\G £ {x,y) - G {x,y)\ < -j— for any x, yen, (3.8) 

\x — y\ a 1 

where C depends only on d, m, \i, fl as well as A and r (if m > 1). 
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Proof. Under the assumptions on A and Q, the estimates \G £ (x,y)\ < C \x — y\ 2 d and 
\V x G (x,y)\ < C\x — y| 1-d hold for any x,y G Q and e > 0. We now fix x ,y G Q and 
r = \x -y \/8. Let fe C™{D{y ,r)), 

u E (x)= / G £ (x,y)f(y)dy and m (x) = / G (x,y)f(y)dy. 
Jn Jn 

Then C £ {u £ ) = £o(uo) = f in Q and w £ = Mo = on <9fi. Also, note that since Q is C 1 ' 1 , 
{ ||V 2 no|| L p(n) < Cp||/|| X p( D(2/0)r)) for any 1 < p < oo, 



[ ||Vuo|U<»(n) < Cp r p \\f\\LP(D(y ,r)) for any p > d 



(3.9) 



(see e. g. [16]). 
Next, let 

w e = u £ - u - eXj (x/e) ^ = £ ( x ) + ^(x), 

where 6 £ G H^{Q) and £ e (fl e ) = £ £ {w e ) in fi. Observe that by (ESj) with V^f e = Pf (x) + 
exj{x/e) and Remark EZQ ||0 e ||#jL(n) < C e\\ V 2 u \\ L 2 [n) < C e\\f\\ L 2 {D{y(hro)) . By Holder's and 
Sobolev's inequalities, this implies that 

\\0e\\v{D{x o ,r)) < Cr\\6 e \\ Lm < C er\\f\\ L 2 (D{yQ , r)) < C er 1+ ^l\\f\\ LP{D{yo , r)h (3.10) 
where q = and p > d. Also, note that since C e (z £ ) = in Q and z £ = w £ on <9f2, 

||-Ze|U»(n) < C\\z £ \\ L ^ {dn) < C e\\Vu \\ L oo( m y (3.11) 
In view of (I3.9I) - (I3.11I) . we obtain 

\\u £ — Uo\\L 2 (D(x ,r)) < \\9e\\L 2 {D(x ,r)) + \\z £ \\L 2 (D{x ,r)) + C ST 2 || Vm || L°°(Q) 

< Pe\W{D{x ,r)) + Cer* || Vito|U»(n) (3-12) 

< C£r 1+ 5-|||/|| iP(D(y0ir)) , 

where p > d. This, together with Lemma [3.21 and (13.91) . gives 

|u e (x ) — tto(a?o) | < Cp 11/11 LP(D(y ,r))- 

It then follows by duality that 

f /" ' Y IP ' i-* 

<^ / |G £ (x ,y) - C (x , 2/)| p ' <iy f <C p er 1_ p for any p > d. 

lJD(y ,r) J 

Finally, since £* (G e (xo, •)) = ^0(^0(^0,')) = in D(y ,r), we may invoke Lemma f3T2l 
again to conclude that 

C f 

\G £ {x ,y ) - G {x ,y )\ < — / \G £ (x ,y) - G {x ,y)\dy 



r JD(y ,r) 

+ C e\\ V y G (x , •)\\L° (D(yo,r)) 

+ Cper^p II VjG o (ar o ,0lli>(2J(j«, ) ro)) 
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where we also used 



l/p 

\V 2 y G (x ,y)\ p dy }> < C p r' 2 \\Go(x , •)IU°°(D(vo,ar)) < Cp r " 



'D(y ,r) 

obtained by the boundary W 2,p estimates on C 1 ' 1 domains [16j. This completes the proof. □ 

As a corollary of estimate (I3.8p . we obtain an 0(e) estimate for ||it £ — Uo||ip(n) for any 
p > 1. In particular, we recover the estimate \\u £ — «o|U 2 (n) < C^||-^||l 2 (q), proved in [17] 
for scalar equations with bounded measurable coefficients satisfying fll.2p - fll.3p . Also see [20] 
for estimates of \\u £ — %||L 2 (n) on Lipschitz domains. 

Theorem 3.4. Suppose that A(y) and Q satsify the same conditions as in Theorem \3.3[ For 
F G L 2 (Q) and e > 0, let u £ G (Q) be the solution of C £ (u e ) = F in Q. Then the estimate 

IK - Wo||w(fi) < Cs\\F\\ LP ( n ) (3.13) 
holds if 1 < p < d and \ = ^ — {\, or p > d and q = oo. Moreover, 

\\u £ - uo|U-»(n) < Ce [In (e~ l M + 2)] x ~* ||F|| L<i(n) , (3.14) 

where M = diam(fi). 

Proof. It follows from the Green function representation and Theorem 13.31 that 

f \ F (y)\ 

\u £ (x) — uq(x)\ < C e j- dy, for any xeH. 

Jn F — 2/1 

This leads to (I3.13P for 1 < p < ti and i = ^ — ^ by the well known estimates for fractional 
integrals. The case of p > d and q = oo follows directly from Holder's inequality. To see 
(I3.14p . we bound \G £ (x, y) — G (x, y)\ by C\x — y\ 2 ~ d if \x — y\ < e, and by Ce\x — if 
\x — y\ > s. By Holder's inequality this gives 

\u £ (x) - u (x)\ < c f JEM_ dy + C sf j^Ldy 

J D(x,e) F ~~ V I Jn\D(x,e) F — i/ 1 

<C£||F|| id(n) + C'£[ln (e- 1 M + 2)] 1 ^||F|| Ld(n) 
<C £ [ln ( £ - 1 M + 2)] 1 ^||F|| Ld(Q) , 
which completes the proof. □ 



3.2 Lipschitz estimates 

In this subsection we give the proof of ( 11. 8ft . As a corollary of ( 11. 8ft . we also obtain an 0(e) 
estimate for u £ - u - {$ £J - - Pj}^ in Wq P (Q) for any 1 < p < oo. 

Recall that D(r) = D(x ,r) = B(x ,r) fl Q and A r = A(x ,r) = B(x ,r) PI dfl, where 
x G f2 and < r < r . Throughout this subsection we will assume that Q is a bounded 
C 2 '^ domain for some r\ G (0, 1) and A = A(y) satisfies conditions ( 11.21) . (II. 3ft and (11.41) . 
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Lemma 3.5. Suppose that u £ G H l (D ir ), u G C 2,p (D ir ) and C £ {u £ ) = £o(u ) in D ir , where 
< p < r\. Also assume that u e = Uq on A 4r . Then, if < e < r, 



<~5^T / \u £ -uo\dx + Cer~ x \\Vuo\\ L <x, {Dir) 



~~ — \ ' ~\ ' ~~~ 1 1 L 00 ( D ~) 

dxi dxi 1 £,J ' dxj r 

j \u £ - w 1 dx + Cer' 1 \\ Vu \\ L 
+ Ce In [e-V + 2] || V 2 u \\ L »° { D ir ) + Cer p \\ V 2 u \\ 
where $ e = denotes the matrix of Dirichlet correctors for C £ in Q. 

Proof. We begin by choosing a C 2,v domain D such that D^ r G D G D^ r . Let 



/8 _ p/6f\ . ^"O 

dx 



3 



Note that w £ = on A 4r . Write u> £ = itfg + wi in D, where We G H^(D) and ^(we 1 ') = 
£ e (u> e ) in .D. Since C £ {wf^) = in _D 3r and = «j e = on A 3r , it follows from the 
boundary Lipschitz estimate in [21 Lemma 20] that 



liv<4 2) iu~ ( zv)<-^i / kfl^ 



g 



£>2r 



/ \w £ \ dx + Cr 1 \\wW\\ L c a ( D2r ) 

JD 2 r 
C f 

-~d+T / \u £ -UQ\dx + Cer- l \\Vu \\ L ^ {D2r) + Cr~ 1 \\w^ ) \\ L ^ {D2r) , 
r Jd 2t 

where we have used the estimate ■ — Pf ||oo < Cs in Proposition 12.41 This implies that 



C f 

||Viu B ||i,ocp r ) < -j-t / |M e -M Ma; + C^ _1 ||VMo|| L oo (D2r) + C||Vw( 1) || L o 



where we have used ||u>! Hl 00 ^,.) — C^UVwe ||x°°(D2r)- Thus, 

\®e,jl ■ —\\L~iDr) 



(D 2r ), 



dxi dxi 1 6,3 i dxj 

- Zd+x [ \ue-u \dx + Cer~ l \\Vu \\L°°(D 2r ) ( ' > ' 10) 



r" 



D 2r 

+ Ce\\V 2 u \\ L ~ {D2r) + C\\VwP\\ L ~ a ^ )■ 



To estimate VtOe on D 2r , we use the Green function representation 

wP(x)= ld e (x,y)C e (w e )(y)dy, 
Jd 
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where G £ (x,y) is the matrix of Green functions for C £ in the C 2,v domain D. Let 

fi(x) = eF kij (x/e) ® ^° + a {j (x/e) [$ 6>fc - P k ] ■ 



OXjOXk OXjOXk 

where we have suppressed the superscripts for notational simplicity. In view of ( 12. 8p . we 
obtain 

wP(x) =- J JL{G g (x,y)} • {My) - /,(*)} dy 

+ J G e (x, y)a ij (y/e) — [$ 6 , fc - P h - e\k {y/e) ] ■ ^ ^ dy. 

It follows that 

\VwP(x)\< [ \V x V y G £ (x,y)\\f(y)-f(x)\dy 

+ C\\V 2 u \\l~(d, t ) [ \V x G £ (x,y)\\Vy[$ £ j - Pj - E Xj (y/e)]\dy. 
Jd 



(3.17) 



To handle the first term in the right hand side of A3 . 1 Tj) . we use \V x V y G £ (x, y)\ < C\x — y\ d 
and the observation that 

\\fh^(D 4r ) < Ce\\V 2 u \\ L c iD4r) , 
\f(x) - f(y)\ <C\x-y\" {e l - p \\V 2 u \\ LOO(Dir) + e\\ V 2 u Q \\ c ^ { D ir )} • 

This yields that 

jV x V„G e (z, y)\\f(y) -f(x)\dy 

dy 



D 



°° (Air) /_ 

JD\B(x,s) 



x — y\ d 

dy 



+ C {e 1 -^ V 2 M || L - (D4r ) + e\\V 2 u \\c°, P{Dir )} [ 

JDnB(x,s) 

<C£ln[£- 1 r + 2]||V 2 Wo|U-(^ r ) + C£ 1 +^||V 2 Mo||co,P(z? 4r )- 



\x - y\ d ~ p 



Finally, using the estimates \V x G £ (x,y)\ < Cdist(y, dfl)\x — y\~ d and \V x G £ (x,y)\ < 
C\x — y\ l ~ d as well as estimates in Proposition I2.4[ we see that the second term in the 
right-hand side of (13.171) is bounded by 

I JD\B(x,s) \ x — y\ J DnB(x,e) F ~ V\ J 

<Ce\n[s- 1 r + 2]\\V 2 u \\ L o O{D4r) . 
Thus we have proved that 

HV^IUoo^) < C^^ln^-V + 2]||V 2 w |U-(^) + Ce^^HV^oll^PC^,-)- 
This, together with (13. 16ft . completes the proof of ( 13. 15ft . □ 
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We are now ready to give the proof of the estimate (11. 8p . 

Theorem 3.6. Suppose that A(y) satisfies conditions U.S\) . U.3\) and fli.^| ). Let Q be a 

bounded C 2,r? domain for some r\ e (0, 1). Then 

\1\G7(*,V)} - ■ £{<#V.v)}| S C£ ' n| ; 1 ^ l + 21 . (3-18) 

for any x,y G Q, where C depends only on d, m, fi, X, r and Q. 

Proof. Fix x , y E Q and r = \x — y \/8. We may assume that e < r, since the case e > r 
is trivial and follows from the size estimate of \V x G E (x, y)\, \V x Gq(x, y)\ and ( 12. 121) . 



Let u e {x) = G £ (x,y ) and u (x) = G (x,y ). Then C £ (u £ ) = C (u ) = in D 



D(xo,4r) and u e = uq = on A^. = A(xo,4r). Note that by Theorem 13.31 we have \\u £ — 
«o|U°°(£> 4r ) < Cer l ~ d . Also, since Q is C 2 '' 7 , we have || Vu ||L°°(D 4r ) < Cr l ~ d , \\ V 2 ii |U°°(£> 4l .) < 
Q r -d anc [ || V 2 Uo||c°p(D 4r ) < Cr~ d ~ p . It then follows from Lemma [3751 that 

»S-!-f*"}-gll-( D ^^H,-V + 2]. 

This finishes the proof. □ 

As a corollary of the estimate (13.181) . we obtain an 0(e) estimate (up to a logarithmic 
factor 'dp ^2) for u £ - u - - P^fg in W^ p (n). 

Theorem 3.7. Assume that A(y) and Q satisfy the same assumptions as in Theorem \3.b\ 
Let 1 < p < oo. For F e L p (n) and e>0, let u £ e W£*(Q;) and C £ (u £ ) =F into. Then 



\u. 



«o - {<; - P}}^\\wt»W < C M He- l M + 2]}^\\F\\ LP{n) , (3.19) 



where M = diam(fi) and C p depends only on d, m, p, fi, X, r and Q. 
Proof. We will show that for any 1 < p < oo, 

11 tr " ^ { *« } ' tr iu " n > - Ce {lll|£ " M + 2|}4 ' H ' l|F|U " n »- (X20) 

This, together with (12. 12j) and the estimate || V 2 w ||LP(n) < C||P|Up(n) for 1 < p < oo, gives 
(I3.19p . To see ( )3.2Up . we use Theorem 13.61 as well as estimates on \V x G £ (x,y)\ and |V$ E | to 
deduce that 

, dut d '■ 1 " ' 



WJ}-^|<^ / tf«(*,i/)l/(v)l<&, 

3 J Q 



where 

*.(»,»)=< r |l- rr^ [e " 1|l " l,l+21, i !! I " !, !- £ ' (3.2i: 

\x — y\ , it \x — y\ < e. 
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Note that 



sup / K £ (x, y) dy + sup / K £ (x,y)dx<Ce{la[e- 1 M + 2)y 
ie!) Jn y&n Jn 



This gives (I3.20p in the case p = 1 or oo. Thus, by the M. Riesz interpolation theorem, it 
suffices to prove the estimate for p = 2. 

Let w £ = u £ — u — — Pj}^r~ an d S(x) = dist (x, dQ) . We may deduce from 

Propositions O and El that 



/ \Vw £ \ 2 dx < Ce / \V 2 u \\Vw £ \dx + Ce / [5(x)] _1 |V 2 u | \w e \ dx 
Jn Jn Jn 

< Ce\\V 2 u \\ L 2 {n) \\Vw £ \\ L 2 {n) , 



(3.22) 



where, for the last inequality, we have used Holder's inequality as well as Hardy's inequality 
||[5(x)] _1 w e ||i2(Q) < C\\ Vw £ ||z,2(q). The desired estimate follows easily from (13. 22 p . □ 

Let G* £ (x,y) = (G* a/3 (x, y)) denote the matrix of Green's functions for £*, the adjoint of 
C £ . By Theorem 13.61 

" ^WW) ' £;{Gf{*.V)}\ S CE, 'X l * y } yl + 2] - (3-23) 

where $* denotes the matrix of Dirichlet correctors for C* in Q. Since G* a/3 (x, y) = G^ a (y, x), 
we obtain 

l>.)l-^i.^(,,,)H < C£ll ' l ;'ly l + 21 (3.24) 

This leads to an asymptotic expansion of the Poisson kernel for £ e on f2. 
Let {h a P(y)) denote the inverse matrix of (ni(y)rij(y)af-) . 

Theorem 3.8. Suppose that A(y) satisfies conditions ( UUP , ( f and ( fi.^| j. Z-et P £ (x,y) = 
(P £ Q/3 (x,y)) &e i/ie Poisson kernel for C £ on a C 2,v domain Q. Then 

Pf(x,y)=P^(x,y)uf(y) + Rf(x,y), (3.25) 

where 

"f(y) = h"(y) ■ g^{K7(y)} ■ Mv) ■ THMnsMagiy/e), (3.26) 



\Rf(x,y)\ < l - — 1 — ,/' J for any x efl and y e dQ, (3.27) 

\x — y\ d 

and C depends only on d, m, fi, X, r and Q. 
Proof. Note that 



Pf(x,y) = ^{v)a${vle)^{G?{x, y )} 
d 

{G^(x,y)} ■n i {y)n j {y)alHy/e), 



(3.28) 



dn{y) 
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since G e (x, •) = on dQ. By (I3.24p . we obtain 



Celnle- 1 ] 


x - 


-y\ 


+ 2] 




x - 


-y\ 


d 



(3.29) 



Celn^- 1 ! 


x - 


-y\ 


+ 2] 




x - 


-y\ 


d 



In view of (13. 28ft (with e = 0), we have 

P?(x,y)h*{v) = -^{Gr(x,v)}- 
This, together with (I3.29p . gives 

\Pf(x,y)-P^(x,y)uf(y)\< 

for any x G VL and y G dQ, where oo £ (y) is defined by (I3.26p . □ 

With the asymptotic expansion for the Poisson kernels at our disposal, we may approx- 
imate the solution of the LP Dirichlet problem: C £ (u £ ) = in Q and u £ = f £ on dfl by 
the solution of the homogenized system with boundary data u> £ f £ . As mentioned in the 
Introduction, the case f £ {x) = f(x,x/e) with f(x,y) periodic in y is of particular interest. 

Theorem 3.9. Assume that A(y) and Q satisfy the same assumptions as in Theorem \3.8i . 
Let C £ {u £ ) = in Q and u e = f £ on dQ. Then for any 1 < p < oo, 

(■ r> -> 1/p 

\\u £ -v £ \\ LPm <C{e(ln[e- l M + 2]) ] \\f £ \\ LP{d n), (3.30) 

where C (v £ ) = in Q and v £ = u £ f £ on dQ, with u £ defined by A3.26\) . 
Proof. Since 

u £ (x) = P £ (x,y)f £ (y)dy and v e (x) = P (x,y)u E (y)f £ (y) dy, 
Jan Jan 

it follows from Theorem 13.81 that 

\u e {x) -v e (x)\ < / \R £ (x,y)\\f e (y)\dy. 
Jan 

Using 

\Re(x, y) I < C{\V y G £ (x, y) | + \V y G (x, y) |} < C5(x)\x - y\-\ 

where 8{x) = dist(x,<9fi), we see that f gn \R £ (x,y) \ dy < C for any x G Q. It then follows 
by Holder's inequality that 

\u £ (x)-v £ (x)\ p <C f \R £ {x,y)\\f £ {y)\vdy. (3.31) 
Jan 
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Now, by Theorem 13.81 as well as the estimate \R £ (x,y)\ < C\x — y\ l d , we obtain 



/ 

Jn 



dx f ... | | _ 

R £ {x,y)\dx < C I -j^ + Ce ^ — *J ' ~ J dx 

QnB(y,e) \ X ~ V\ Jd\B(y,e) 



lnfe" 1 


\x-y 


1 + 2] 




x — y\ 


i 



< C£{ln[e- 1 M + 2]} 2 , 

for any y G dQ. This, together with (I3.3ip and Fubini's Theorem, gives ( 13.30!) . □ 
We end this section with another approximation result. Note that by ( 13.201) . 

\\V(C £ )-\F) - V$ £ • VOCor 1 ^)!!^) < Ce{ ln^M + 2]} 4 '^' \\F\\ LP{n) , (3.32) 
for 1 < p < oo. By duality this gives the following. 

Theorem 3.10. Assume that A(y) and Q satisfy the same assumptions as in Theorem \3.8\ . 
For f = (/f) G L 2 {VL) and e > 0, let u £ G H%(£1) and C £ (u £ ) = div(f) in Vt. Then if 
f G L P (Q) for some 1 < p < oo, 

\K - v £ \\ LP(n) <Ce{ He^M + 2]} 4 '^' ||/|| LP(n) , (3.33) 
where v £ G Hq(Q) and C (v £ ) = div(F e ) with 



F ^ x) = f ^ x) ^{ 



*/3a 



dXj 



£.1 



3.3 An asymptotic expansion for V x V y G £ (x, y) and its applications 

In this subsection we derive an asymptotic expansion for V x V y G £ (x, y). As its applications 
we obtain asymptotic expansions for (d/dxi)(£ £ ) (d/dxj) and the Dirichlet-to- Neumann 
map associated with C £ . 



Theorem 3.11. Suppose that A(y) satisfies conditions lil.fy) , M.3}) and (fi. J$ . Let Vt be a 

C 3,v domain for some rj > 0. Then 



dxidyj £ dxi £ ' dx k dyi dy 



(3.34) 

^ i ■;■ id : -\.r — // l\ 



Celn [e- 1 


X 


-y 


+ 2] 




x — 


y 


d+l 



for any x, y G Q, where C depends only on d, m, n, X, r and fl. 

Proof. Fix x ,y G Q. Let r = \x — yo\/S. Since \V x V y G £ (x, y)\ < C\x — y\~ d , it suffices to 
consider the case e < r. Fix 1 < /3 < m and 1 < j < d, let 

dG af} 
dyj 

f) f)C acr 

uaM = W^ ]iyo) '^yT {x,yo) 
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in D 4r = D(xq,At). In view of (I3.24[) we obtain 

\\u £ -u \\ Lo . {Dir) <Cer- d \n[e- l r + 2\. (3.35) 
Since fl is C 3,v , we have || Vwo||l°°(d 4i .) < Cr~ d , 

|| W)|U~(iv) < Cr- d ~ x and || VV)||c°-(D 4r) < Cr^- 1 ^. (3.36) 
By Lemma 13.51 estimates ( 13.351) and ( 13.361) imply that 

||^-A/^n < Ceh^r + 2\ 

11 dxi dxi X £ ' k * dx k Uaa{Dr) ~ r d + l 

This gives the desired estimate (I3.34D . □ 



Let u £ G H^(Q) and C £ (u £ ) = div(/) in ft, where / = (//) G L 2 (ft). Then 



9< 
dxi 

£fl wj ), wh ere 

= p.v. £ ^-{Gf{x, y)}g{y) dy. 

It is known that if Q is C 1,r) and A(y) satisfies conditions (11 .2p -( |TT4l) . operators T £ fj are 
uniformly bounded on L P (Q) for 1 < p < oo, and of weak type (1, 1) [5]. 

Theorem 3.12. Suppose that Q and A{y) satisfy the same conditions as in Theorem \3 . 1 1[ 
Let 1 < p < oo and g G L P (Q). Then as e — > 0, 



(3.37) 



^ in L 9 (ft), 

/or any 1 < q < p. 

Proof. Let S e (g) denote the left hand side of (13.37)) . By uniform boundedness of ||T"^-|| 
||V$ e ||oo and ||V$*||oo, we see that the operators S £ : L P (Q) — > L q (Q) are uniformly bounded 
for 1 < q < p < oo. As a result we may assume that g G C 1 (IR d ). 
Note that S £ (l) = and 

= *fa - *(.))<»> = p-v. / *.(*,»){,(») - rf,)} *, 

JO 

where, by Theorem 13.111 and the estimate | W x W y G e (x, y)| < C|x — the integral kernel 
K £ (x, y) satisfies 



1 


elnfe- 1 ! 


x - 


"2/1 


+ 2] 


-y\ 




x - 


-J/l 





|-Ke(z,y)| < C'min 
It follows that if e < (1/2), 

n , u ^ „„ „ /" dy „ ., /" hit -1 \x — V I + 2] , 

p E S W < C </ c i (!2) / -I rj-r + Ce \\g £00(n) / : ^— 

< Cy/e\ lne|||p||cn(n)- 

Hence, || ^(c/) || z,«?(o) — > as s — > 0. This completes the proof. □ 
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Finally, we consider the Dirichlet-to-Neumann map A e associated with the operator C £ . 
Let / G H 1 / 2 (dQ,) and u £ G H 1 ^) be the solution of C £ (u £ ) = in and u £ = f on 
dtt. The map A £ : H l l 2 {dtt) -> H- l / 2 {dtt) is defined by A £ (/) = g£. It is known that 
A £ : W 1 ' p (dQ) — )■ L p (dVt) is uniformly bounded for 1 < p < oo, if f2 is C 1 ^ [21]. In the case 
that O is Lipschitz, the map is uniformly bounded for l<p<2 + <5ifm = l [22], and for p 
close to 2 if m > 1 [23J. For simplicity we will assume m = 1 and A* = A in the rest of this 
subsection. 

Let 

d 

u E (x) = ni(x)nj(x)aij(x/e) ■ [n h {x)m{x)a M ] ■ — {$ £>s } ■ n s (x). 
Then ||w £ ||loo(^) < C. It follows from Theorem 13.111 that 



9 9 

{Pe(x,y)}-u} e (x) ■ {P (x,y)} ■ U) £ 



1 du £ {x) du (x 
< 



Celn^ 1 ! 


X 


-y\ 


+ 2] 




x — 


y\ 


d+l 



(3.38) 



for any x, y G dQ. Using $ £j fc(x) = Xk on <9f2, one may show that ninjdijU £ (x) = rikA £ (xk). 

Theorem 3.13. Suppose that Q and A(y) satisfy the same conditions as in Theorem [X721 
We further assume that m = 1 and A* = A. Let f G if^dfi). Then 

df 

Ae(/) - rh-^-Aeixj) + uj £ [/AoK) - Ao(w e /)] 

(3.39) 

+ LO £ ni— [A (u s Xj) - XjA (u £ )] 0, 
m L q (dQ) for any 1 < q < 2, where J^- = (^igfr — n j~£r)f- 

Proof. By a linear change of variables Theorems 15.11 and 15.21 in Section 5 continue to hold 
for A = A . It follows that ||A (o; e /) - fA (u £ )\\ L 2 {m) < C||/|| H i( 9 n) and ||A (o;^-) - 
XjA (u £ )\\LP(gn) < C p \\uj £ \\lp^qq) < C p for any 1 < p < oo. Hence, by Holder's inequality, the 
left hand side of (13.391) . as an operator, is uniformly bounded from ff 1 (9fi) to L q (dQ) for 
any 1 < q < 2. Consequently, it suffices to show that for / G C 2 (dQ), the left hand side of 
(I3.39P goes to zero in L°°(<9f2), as e — > 0. 

To this end, recall that A e (/) = where u £ (z) = J n P £ (z,y)f(y) da(y) for z G O. 
Write 

u £ (z)-f(x)= [ P £ (z,y){f(y)-f(x)}da(y) 
Jan 

Pe(z,y) \f(y) - f(x) -Ui{x) d{ . (yj-Xj) \ da{y) 
an I otij{x) j 

+ / P e (z, y)ni(x)j-^ ■ - arj) da(y) (3.40) 
./an otij{x) 

Pe(z,y) \f(y) - /(x) -7ij(x) ^{ • (yj-Xj) \ dcr(y) 
an I otij{x) j 
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Since \ f(y) — f(x) — n^x) ■ df /dt^^x) ■ (yj — Xj)\ < Cf\y — x\ 2 for x, y G dfl, it follows by 
taking derivatives in z and then letting z — > x in (I3.4(jp that 

f d f df 1 

A E (f)(x) = J ^ {P E (x,y)} U(y) - /(^-"i^ ^Tj " j dcr (v) 

df 

ot ij (x) 

In view of (13.38}) as well as the estimate \V x P £ (x,y)\ < C\V x V y G £ (x, y)\ < C\x — y\~ d , we 
obtain 

df 

A E (f)(x) = I E (x) +n i -^-A E (x j ) 

f d ( df 1 

+ / UJ e(x)——-{PQ(x 1 y)}u £ (y) If (y) - f(x) -ni(x)-—-^ ■ - xj) \ da(y) 
Jon ovq{x) \ dUjix) J 

df df 
= I E (x) + ni—A E (xj) + u £ [A (u f) - fA (uj E )] + ^^gf- [x 3 -A (w 6 ) - A (u; £ x j )] , 

where the term I £ (x) satisfies 



ln[e- 


-i 


x — 


y 


+ 2] 




y 


— X 


d- 


i 



\h{x)\<C f le / -V ^^d (7 (y) + 

I Jdn\B(x,e) 

< C / £[ln(e- 1 M + 2)] 



Id— 2 



'B(a:,£)n<9n 12/ ~~ ^1 

This gives the desired estimate. □ 



4 Asymptotic behavior of Neumann functions 

Throughout this section we will assume that A(y) satisfies conditions (ll.2p . (II. 3p and (II. 4p . 
Under these conditions one may construct a matrix of Neumann functions N £ (x,y) = 
(N°P{x,y)} in a bounded Lipschitz domain fl such that 

C £ {N £ ^;y)} = e^ y (x) in fi, 
< ^{JVf^jH-elffll- 1 on90, (41) 

/ N?(x,y)da(x) = 0, 
^ Jan 

where JVf (x, y) = {N^{x, y),..., N^(x, y)) and = (0, . . . , 1, . . . , 0) with 1 in the (3 th 
position. Let u £ G H l (Vl) be a solution to C £ (u £ ) = F in f2 and |^ = g on <9fi. Then 

Ue ^ ~WT\J dn U£ = J Q N ^ Xj y)F ^ y) dy + J m N ^ Xj y)9i ~ y) da ^' (4 ' 2) 
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Under the additional assumption that A* = A and Q is C 1,v for some 77 G (0,1), it was 
proved in [2T] that 

\N s {x,y)\< 



\V x N{x,y)\ + \V y N £ {x,y)\ < 
\V x V y N E (x,y)\ < 



c 


\x 


-y 


d-2 




c 




\x 


-y 


d-l 




c 




\x 


- y 


d 



(4.3) 



for any x, y G ft. The goal of this section is to establish the asymptotic estimates of N £ (x, y) 
and V x N £ (x,y) in Theorem II .21 



4.1 L°° estimates 



The goal of this subsection is to prove the estimate (11.91) . We also obtain an 0(e) estimate 
(up to a logarithmic factor) for ||ii e — wo||z>(f2) for solutions with Neumann conditions. Recall 
that D r = D r (x , r) = B(x , r) Dfl and A r = A(x , r) = B(x , r) fl dfl for some x G fl and 
< r < r . We begin with an L°° estimate for local solutions. 

Lemma 4.1. Let fl be a bounded C 1,n domain for some rj G (0,1). Let u e G H l (Dz r ) and 
uq G W 2 ' p (D Sr ) for some p > d. Suppose that C £ [u £ ) = C (u ) in D 3r and |^ = |^ on A 3r . 
Then, if < e < (r/2), 



\U £ — Uo\\L°°(D r ) 



C f 

- ~d / lus-uol + Celnls^r + 2}\\Vu \\ L ™( D:ir) 
r JD 3r 

+ C p er l ~v\\y 2 Uo \\ LP{D , ir) . 



(4.4) 



Proof. By rescaling we may assume r = 1. Choose a C 1,v domain D such that D 2 C D C .D3. 
Let 

w £ = u £ (x) - u Q (x) - EXj(x/e)-^-. 
Using Proposition 12.21 with vfj = Pj(x) + eXj(x/e), we see that 



(4.5) 



where b^ k (y) = F^ k (y) + '{y)Xk \v) * s a bounded periodic function. Also, by a direct 
computation, we have 



dw £ \ a ( du £ 



dv F 



e ( d 
+ 2 n ^ 



-embff k (x/e)- " 



(4.6) 



dxjdxk 
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(see j20j Lemma 5.1]). Let w £ = wi + We, where 

(«?>(*))• = -e I £-{S?{x,y)} ■ bfl(y/ S ) ■ gjL* (4.7) 



A' 



and N e (x,y) denotes the matrix of Neumann functions for C £ in D. Since |V y A^ e (x,?/)| < 
C|s — y\ 1 ~ d , it follows from Holder's inequality that 



w 



Wl 



l-(d 2 ) < Cp£||V 2 Mo|| LP ( D3 ) for any p > d. (4.8) 



To estimate wf\ we observe that £ e (u4 2 ' 1 ) = in D and 



on <9-D. Let Wg 2 ' ) = -u4 21 ^ + u>£ 22 \ where 

For each x E D, choose x G <9.D such that |sc — x| = dist(x,dD). Note that for y G <9-D, 
|y — x\ < \y — x\ + \x — x\ < 2\y — x\. Thus, \V y N £ (x, y)\ < C\y — x\ l ~ d and 

lM 21, )»i = f I / sg -«.|t) {^(x,„)} ■ {/?(„) -/?(*» <M!/)| 
< C£ /Mi, (s)i 

where /(y) = ($(y)) = (F^fo/e^y)). Since ||/|| L - (A|) < C|| Vu \\l~(d 3 ) and 



I/O) - /(«)| < Ce l \y - x\\\Vuq\\ l ^ {Da) + |y - Vm ||c°.o(d 3 ), 



where < p < r] and we have used the fact || Fj^Hc^Y) < C, it follows that 



H 21 \x)\ < CeHVuolU-^a) / \x-y\ l - d da(y) 

JdD\B(x,e) 



+ C||Vuo|Uco (Z) 3) / |y - x\ 2 ~ d d(j{y) 

J B(x,e)OdD 

+ Ce\\Vu \\ c o, P{D3) [ \y - x^ da(y) 

JB{x,s)ndD 

< Ce\n[e- 1 + 2}\\Vu \\ L ^ (D3) + Ce 1+p \\Vu \\ c0 , P{D , i) . 



By Sobolev imbedding, this implies that 



l«4 21) IU~p 2 ) < Celnie- 1 + 2]\\Vu \\ L o a{D3) + C p e\\V' 2 u \\ LP{D3) (4.11) 
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for any p > d. 

Finally, to estimate wf 2 \ we note that C £ {w {22) ) = in D 2 and -^{wf^} = §j£-fjj = 
on A 2 . It follows from [2TJ Theorem 3.1] that 

\\w^\\ T ^nA<C / \w^\dx 



J £ '||L<»(Oi) ^ ^ / 

'£>2 

< C / \u £ - U \ dx + C£||Vtio||L°°(D 2 ) 

+ CH^Wllioc^) + C\\w {2l) ||l-(d 2 )- 

Hence, 

||«e - «o|U°°(.Di) <C / | M e - U \ dx + || Vu || L°°(D 2 ) 
-/D 2 

+ cWwPWlo*^ + c||w( 21) || LK > (Z)2) . 

This, together with (g3J and fHTTjl . gives (Ojl . □ 

The next lemma on the traces of fractional integrals is known. We provide a proof for 
the sake of completness. 

Lemma 4.2. Let Q be a bounded Lipschitz domain. Let g G Cq(Q D B(y ,r)) for some 
y E Q and 

g(y) dy 



u\x-y 



d-l 



u(x) 

Then \\u\\ L 2 {m) < Cr 1/2 |MU2 (n) . 

Proof. By the well known estimates for fractional and singular integrals, ||Vw||i2(n) < 
C||<?IU 2 (n) an d IMU 2 (n) < C||fl'||i>(n), where p = j^. Choose a vector field v G C^°(M d , R d ) 
such that <w,n>>co>0on dQ and \v\ < 1, \Vv\ < C/r in ~R d . It follows from integration 
by parts and Holder's inequality that 

cq \u\ 2 do < I < v,n > \u\ 2 do < Cr^ 1 \u\ 2 dx + C / |tt||Vu|rfx 
Jan Jan Jq Jq 

< Or' 1 I \u\ 2 dx + Cr I \Vu\ 2 dx 

Jq Jq 

< Cr^WgWl^ + Cr\\g\\ 2 L2{n) 

where we have used the fact supp(g) C B(y ,r). This completes the proof. □ 
We are now ready to prove the estimate (11.91) . 

Theorem 4.3. Suppose that A(y) satisfies conditions U.fy) . M.'J\) . Ijl.4\ ) an d A* = A. Let 

Q be a bounded C 1 ' domain. Then 



\N £ (x,y)-N (x,y)\ < (4.12) 

\x — y\ a 1 

for any x, y G ft, where C depends only on d, m, n, X, r and Q. 



Celn^" 1 


X 


-y 


+ 2] 




x — 


y 


d-i 
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Proof. By rescaling we may assume that diam(fi) = 1. Fix xq,Uq G Q and let r = |x — 2/o|/8- 
Since \N £ (x ,y )\ < Cr 2 ~ d , we may assume that e < r. For g e C^(D(y ,r)) and e > 0, let 
= /n iV ^ x ' J/Mf) Tnen £ e (u e ) = 5 in fj£ = -pjf / n 9 on <9fi and J 9n u £ = 0. It 

follows that C £ (u £ ) = C (u ) in and ^ = |^ on <9f2. Let u> £ = u £ (x) — u Q (x)— sxj (x/e)^. 
As in the proof of Lemma 14. 1[ 

and 

'dw £ \ a _e ( d d \ / x^uol j,afl/ , s d 2 u$ 



dXjdXfr 

on <9f2. Now, write u> £ = 9 e + z E + p, where C e {z e ) = C £ {w £ ) in fl, J n z £ = 0, 

and p = t^jt J* an {u> £ — z £ } is a constant. Note that ||V2 £ ||i,2(q) < Ce\\ V 2 Wo||,L 2 (n) < 

C £ IMU 2 (fi)- Since Jo^e = 0, by the Poincare inequality, we obtain ||z e ||£p(n) < Cte||<7||£2(n), 
where p = -r^. It follows by Holder's inequality that 

-—-ii 

\\Za \\tf(D(x ,r)) < Cr 2 p\\z e \\ L p(D(x ,r)) < Cer\\g\\ L 2 {Q) . (4.13) 
Next, to estimate 9 £ , we observe that C £ {6 £ ) = in Q, f dQ £ = and 
86 £ e ( d 



■DM' 



- rii— n 

ov £ 2 \ oxj 

Using a duality argument and L 2 estimates for the Neumann problem in [22J, we may deduce 
that ||^ e ||z,2(af}) < C^|| Vmo||l2( 9C ) (see the proof of Theorem 5.2 in [20]). By the square 
function estimate for the L 2 Dirichlet problem [20], this implies that 

It follows from the Sobolev imbedding that ||0 £ ||,LPi(n) < Ce\\ Vwo||L 2 (,9n), where p% = jzt. 
By Holder's inequality, this gives 

||^||jy»(D(x ,r)) < Cer^\\Vu \\ L 2 {m) . (4.14) 

Since 

\Vu (x)\ <C 



\g{y)\ dy 



Id— 1 ' 



'n F -y\ 

we may invoke Lemma [4.21 to claim that ||Vuo||x, 2 (n) — C?" ||<?||z 2 (fi)- In view of (14. 14ft we 
obtain ||0 £ ||.L 2 (£>(a; o ,r)) < Cer \\ g \\l 2 (q)- This, together with (I4.13P and the observation 



\p\<C {e\Vu \ + \z £ \}da < Ce\\g\\ L ^), 
Jan 
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gives ||t« e ||La(D(* ,r)) < Cer\\g\\ L 2 (n) . It follows that 

1/2 



(4/ |m £ -m | 2 ) < Cer 2 z d \\g\\ L 2 {n) . 

L r JD(x ,r) J 



r JD{x ,r) 

Since || Vwo|U°°(D(xo,r-)) < CVT"||flr||i2(n) and || V 2 u ||iP(n) < C||#|| iP(n) , by Lemma ED we 
obtain 

|u e (x ) - i*o(ar ) I <C p er » hx[e~ r + 2}\\g\\ LP{n) , 
where p > d. By duality this gives 

i- d [ \N E (x ,y) - N (x ,y)f dy) '* < C v er x ~ d lnfe'V + 2]. 

L r JD(y ,r) ) 



Finally, since 



d .{iV £ (a;,y)} = -^— {iV (x,y)} = --^| onffi, 



|Vj;iVo(a;, y)| < C\x — y\ l ~ d and || V 2 A^ (x , y) \\Lp(D( yo ,r)) < Crp~ d , we may invoke Lemma ED 
again to obtain 

C f 

\N £ (x ,y ) - N (x ,y )\ < — / \N e (x Q ,y) - N Q (x Q ,y)\dy 

r JD{y Q ,r) 

+ Cer 1 - d \n[e- 1 r + 2] 
< Cer 1 - d ln[e- 1 r + 2]. 

This completes the proof. □ 

Remark 4.4. It is not clear whether the logarithmic factor in ( 14. 12ft is necessary. Also, in 
view of Theorem 13.31 on Green's functions, it would be interesting to show that the estimate 
( I4.12p holds for scalar equations with no smoothness condition on the coefficients. 

As a corollary of Theorem 14.31 we obtain an 0(e) (up to a logarithmic factor) estimate 
for ||« e - w |U«(n)- 

Theorem 4.5. Suppose that A(y) and Q satisfy the same assumptions as in Theorem \4-3\ 
Let 1 < p < oo. For e > and F e L p (tt) with J n F = 0, let u £ e W l ' p (Q) be the solution 
to the Neumann problem: C £ (u £ ) = F inVl, |^ = and J dQ u £ = 0. Then 

\K - «o|U«(0) < Ce ln^M + 2] \\F\\ LP(n) (4.15) 
holds if 1 < p < d and \ = \ ~ \> or P > d an d q = oo, where M = diam(fi). Moreover, 

\\u E -uo\\L°°(n) < Ce[He~ l M + 2)] 2 ^\\F\\ Ldm . (4.16) 
Proof. Note that by the estimate (14.121) . 



\u £ (x) -Uq(x)\ < I \N e (x,y)-N (x,y)\\F(y)\dy 

<Celn(e"M + 2)/ ». 

Jn F y\ 

The rest of the proof is the same as that of Theorem 13.41 □ 
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4.2 Lipschitz estimates 

In this subsection we give the proof of the estimate (11.101) . We also establish an 0(e*) 
estimate for u £ - u - {\t> e j - -Pjfj 1 in W ljP (Q) for any t E (0, 1) and 1 < p < oo. 

Lemma 4.6. Le£ Q be a bounded C 1,v domain for some rj E (0, 1] and 

f d 2 

u e (x)=p.v. J^^-^-{N e (x,y)}f(y) dy 

for some 1 <i,j < d. Then 

||«e|U-(n) < C{ ln^M + 2] + M"}||/|| i0 o (n) + Ce^H^f), (4.17) 
where M = diam(f2) and H £r) (f) = sup { ^^-y\^ '■ x,y EVt and \x — y\ < s}. 
Proof. For x E Q, choose x E dQ such that \x — x\ = dist(x, dQ). Note that 

u £ {x)= j ^-{N £ (x,y)} ■ {f(y) - f{x)}dy 

f 9 

+ / {nj(y)-n j (x)}-—{N e (x,y)}'f(x)da(y), 

JdU OX{ 

where we have used the fact f m N £ (x,y) da(y) = 0. This, together with the estimates 
\V x N £ (x,y)\ <C\x-y\ 1 - d and \V x V y N £ (x,y)\ <C\x-y\~ d , gives 

\u, ii <r f \f(v)~f( x )\ a„ i riifii / da ^ 

JQ.\B(x,e) \ x ~ V\ J\y-x\<e \V ~ x \ 

< C||/|U» (n) Infe^M + 2] + Ce"H e)1t (f) + C\\f\\ L - (sl) At*. 

This completes the proof. □ 

The following lemma provides a Lipschitz estimate for local solutions with Neumann 
boundary conditions. 

Lemma 4.7. Let Q be a bounded C 2,v domain for some rj E (0,1). Suppose that u £ E 
H l (D 3r ), u E C 2 ' r '(D 3r ) and C £ (u £ ) = C (u ) in D 3r . Also assume that §^ = §^ on A 3r . 
Then, ifO < s < (r/2), 

dxi dxi 1 e,J ' dxj 

^ f l, „, i i /t — -i i„r^-i n/r i oiiiY7„. ii _ ._ (4-18) 



C f 

- ~d+T / \ Ue ~ u °\ dx + C^ -1 ln[e _1 M + 2}\\Vu \\L°°(D 3r ) 



'D 3r 

+ Ce^Vln^M + 2]\\V 2 u4 L o» {D3r) + CerHnf^M + 2] || V 2 u \\ c o, P{D3r) 

for any < p < min(r], r), where (^f ,-) denotes the matrix of Neumann correctors for C £ in 
fl and M = diam(O). 
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Proof. By rescaling and translation we may assume that r = 1 and OeDj. Let 



w £ = u £ (x)-u (x)-{^ J -P J "}-^. 



3 

Choose a C 2,v domain such that D 2 C D C D 3 . We now write 

w E (x) = [ N £ (x,y)C £ (w £ )dy + [ N E (x,y)^- da(y) + — L- / w £ 



D JdD 



9Ve \8D\ JdD 



for x G -D2, where N £ (x,y) denotes the matrix of Neumann functions for C £ in D. In view 
of Propositions 12.21 and 12.31 we have w £ = w £ + w £ + c, where c = J d ^ w £ , 

+ J„N £ (x,y) ■a ij (y/e) — {^ e)k {y) - P k (y) - exk(y / e)} • Q^Q^ d V 



2; 



and 

w ( 2 )(x) = e / N £ (x,y) ■ ni(y)F jik (y/e) ■ da{y) 
JdD oyjdy k 

(we have supressed all superscripts for notational simplicity). 
To estimate wi 2 ^ in D 1; we note that £ e ( , u4 2 ' ) ) = in .D and 

9 f / 211 <9 2 wo £ /" / / \ <9 2m o 

— {1^ I = eniF jik (x £)- — — ^- / niF jik (x £) - da on A 

<9z/ £ 1 J ' dxjdx k \dD\JdD dxjdx k 

where we have used §^ = §^ on A 2 . Since ||^{we 2) }|U°°(A 2 ) < Ce|| V 2 wo||l°°(d 3 ) and 

II — {^ 2) }|| c o,p (A2) < Ce l - p \\V 2 u4 L ^ {D2) + Ce\\V 2 u \\ c ^ ( D 2 ), 
it follows from the boundary Lipschitz estimate in [2T| Theorem 7.1] that 

||Vwf || L0O(Ol) < Ce 1 "! W)|U~p 2 ) + C£||V 2 Mo ||co,p(d 2 ) + C / |wf - c| 

-/D 2 

for any constant c. This leads to 

llVWellioop^ < HVw^Hloo^) + HVw^llioopi) 

<C f |w £ |^ + C||Vw«|| Loo(I?2) + C£ 1 ^||V 2 wo|U-(D3) 
+ Ce||V 2 M || 
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(4.20) 



Since |^ - Pf\ < (7eln[e -1 M + 2] by Proposition E31 we obtain 



5 < 5 J xT/ a/3X ^0|| 



<C / l^-uol + Celn^M + ^UVwolUoo^ + Celn^M + ^HV^ollL-CDa) (4 ' 21) 
+ Ce 1 ""!! V 2 w |U-(o 3 ) + II V 2 u ||c°.p(£>3) + C|| Vu>« || L ~ (D2) . 



It remains to estimate Vu^ on D 2 . The first two integrals in the right hand side of 
(14. 19[) may be handled by applying Lemma [4.61 on D. Indeed, let 

f(x) = -sF jik (x/e) ■ d ^° a ij (x/e){^ £ik (x) - P k (x)}- "" 



OXjOXk OXjOXk 

Note that 11/11^(5) < Ce Infe^M + 2] || V 2 m ||l-(d 3 ) and 

< C£ 1 -H n [£- 1 M + 2]||V 2 «|Uo 0(D3) +C'£l n [£- 1 M + 2]||V 2 «o||co.P(i 33 ). 

It follows by Lemma 14.61 that the first two integrals in the right hand side of (14.191) are 
bounded by 

C7elB[e- 1 Af + 2]{e-"||V 2 «o|U-( iJ 3) + ||V 2 Uo||c°>p(c 3 )}- 
Finally, the third integral in (I4.19P is bounded by 

C||V u \\l~(d 3 ) J - \x-y\ d ~ l V ' ( ^ 

Using that |V w {* e ,*(j/) - P*(y) - ex*(y/e)}| < C* min (l, e[dist(w, <9fi)] _1 ) , one may show 
that the integral in (I4.22p is bounded by C£[ln(£ _1 + 2)] . Thus, we have proved that 

Ik^lU^D.) <C£ 1 -nn[£- 1 M + 2]||V 2 M o||L-(D3)+C , £ln[ £ - 1 M + 2]||V 2 no||co.p(D3)- 

This, together with (14.211) . yields the desired estimate. □ 

We are now in a position to give the proof of estimate (ll.lOp . 

Theorem 4.8. Suppose that A(y) satisfies conditions U.S\) . U.3\) and fli.^p . Also assume 
that A* = A. Let Q be a bounded C 2 ' v domain for some r] G (0, 1). Then 

I A {Nf{x , y)} - A { , 3W} . £ {k ,m } | £ ^-;H ; - : m + 2 | (4 23) 

/or any x,y £ Q and p G (0, 1), where C depends only on d, m, p, X, r, p and Q. 

Proof. Since \V x N e (x,y) \ < C\x — y\ l ~ d and | V\I/f -| < C, we may assume that e < \x — y\ 
and p is small. Fix x ,y G Q, 1 < 7 < d and let r = \x — yo\/8. Let = N^(x, y ) 

and «o(a;) = Nq 7 (x, |/ ). Observe that 



£ e (w £ ) = £ Oo) = in £>(£(), r), 

-l^l" 1 ^ onA(x ,r). 



du £ \ ( du \ .v-.i-j 
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Also, note that ||Vu |U«>(r>(xo,r)) < Cr 1 d , || V^olU^OD^r)) < Cr d and || V 2 M ||c .p(D(x ,r)) < 
Cr~ d ~ p . Furthermore, it follows from Theorem 14.31 that 

\\u £ - MolU°=(D(x ,r)) < Cer^lnfe^r + 2]. 
Thus, by Lemma 14 .7\ we obtain 

This completes the proof. □ 

Remark 4.9. Under the symmetry condition A* = A, one has N^(x, y) = N^ a (y, x). Thus, 
it follows from (14.231) that for any x,y G tt, 

|£ {J V^, s)} - A. j{Kl{y)} | < CE "^ + 2 ' . (4.24) 

Fix /3 and j. Let 

<(*) = {^(*, v)} and = • ^ W^*, y)} . 

Note that C £ (u e ) = C (u ) = in Q \ {y} and |^ = |^ = on dQ. We may use Lemma 
inland estimate (I4.24p to deduce that if Q is C 3 > v for some r? G (0, 1), 

|^{^(^)} - ' ^{^(*.»)} ■ £{*5M>I 

C p e 1 ^ln[ £ - 1 M + 2] 1 ' j 

~~ \x-y\ d+l -p 

for any x,i/GO and p G (0, 1), where C p depends only on d, m, /i, A, r, p and Q. 

As a corollary of Theorem 14. 8[ we obtain an 0{e t ) estimate for any t G (0, 1) in W 1,P (Q) 
for solutions with Neumann boundary conditions. 

Theorem 4.10. Suppose that A(y) and Q satisfy the same conditions as in Theorem \4-S\ 
Let I < p < oo. For e > and F G L p (tt) with f Q F = 0, Zet m £ G iy 1,p (f2) 6e the solution 
of the Neumann problem: C e {u e ) — F inQ, ^f- = and J dn u £ = 0. Then 

IK - u - - if }^lk4, P( n) < Cts'WFW^n) (4.26) 

/or any £ G (0, 1), where C t depends only on d, m, fx, \, r, t, p and Q. 

Proof. Since H^—if |U~ (n) < Celn[e- 1 M+2] and || Vu ||LP(n) + || V^oIU^q) < C?||F||ip ( n), 
in view of Theorem 14.51 it suffices to prove that 
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We will prove that estimate ( I4.27P holds for any 1 < p < oo. To this end, note that 
u £ (x) = f n N e (x, y)F(y) dy, by Theorem g^l 



8 {<,} ■ g| < C e ^H^M + 2] / (4,8) 



dxi dxj 



for any p G (0, 1). Note that if p < 1 - t, then e^ln^M + 2] < CeK Estimate fTOTj) 
follows easily from (I4.28p . □ 



5 Leibniz Rules for the Dirichlet-to-Neumann map 

Let A : H 1 ^ 2 (dVl) — > i/ _1 / 2 (9fi) denote the Dirichlet-to-Neumann map associated with 
C = —A. It is known that the estimate || A(/) |)x,p(5Ji) < Cpll/ll w^ian) holds for 1 < p < oo 
if is C 1 [12], and for 1 < p < 2 + 5 if is Lipschitz pH [2TJ, E]. The goal of this section is to 
prove the Leibniz estimates that were used in subsection 3.3. We remark that our basic line 
of argument extends to the general case C = — div(^4(x)V). Indeed, Theorem 15.21 continues 
to hold for if A(x) satisfies conditions (11.21) . IIVA^ < C and A* = A. Under the same 
conditions on A, Theorem 15. II also holds in the scalar case (m — 1). If m > 1, estimate (15.11) 
holds for C l domains. The details will be given elsewhere. 

Theorem 5.1. Let C = —A and Q be a bounded Lipschitz domain. Then 

\\A(fg) - fA(g)\\v m < C\\f\\ H i {dn) \\g\\ L ~ {dn) (5.1) 
for any f G H^dSl) and g G L°°(dQ). 

Proof. Let w be the solution to the L? Dirichlet problem with data fg; i.e., C(w) = in Q, 
w = fg on dfl and (w)* G L 2 (dQ), where (w)* denotes the nontangential maximal function 
of w. Let u and v denote the solutions of the L 2 Dirichlet problem with boundary data / 
and g respectively. Then 

A(fg) - fA(g) - gA(f) = ^-{w - uv}. 

To estimate d(w — uv)jdv in L 2 (dQ), let h be a function such that C(h) = in Q and 
(h)* G L 2 {dVl). Since w — uv = on d£l and 

du dv d f du 
C(w - uv) = 2—— = 2- 



by the Green's formula, we obtain 



f -?-{w-uv}.h = 2 f^L. v ~-2 f A(f).g-h. (5.2) 
Jan ov Jn oxj dxj J dn 



Thus, by duality, it suffices to prove the following estimate 

du dh 



n ® x j ® x j 



< C\\f\\m(dn)\\g\\L°°(dn)\\h\\ L 2(dn)- (5.3) 
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The trilinear estimate (15. 3 p on Lipschitz domains is a consequence of a bilinear estimate 
due to B. Dahlberg [S] (also see [TB] for related work). Indeed, since A(h) = in Q and 
(h)* G L 2 (dQ), it follows from estimate (1.7) in [S] that the left hand side of (15. 3p is bounded 
by 



C\\h\\ L 2 {m) \ [ \((Vu)v)*\ 2 da + [ \V 2 u\ 2 \v\ 2 5{x)dx + [ | Vu\ 2 \Vv\ 2 5{x) dx\ 
{Jan Jn Jn J 



1/2 

, (5.4) 



where 5(x) = dist(x, <9f2). The first two integals in (15. 4 p are bounded by C* j | / 1 1 ^ i (^a) 1 1 f 1 1 z.°° (oo) • 
This follows from estimates || (Vu)* \\L 2 {an) < C \\f ||#i (an) an d f n \V 2 u\ 2 5(x) c 
[T9| 17] as well as the maximum principle < Hfi'lU 00 ^)- 



Finally, to handle the last integral in (15 .4p . we use the fact that \Vv (x)\ 2 5(x)dx is a 
Carleson measure whose norm is less than C\\g\\ 2 BMO (dn) P^l - This implies that the integral 
is bounded by CIKVm)*!!^^)!!^!!^^)- Estimate (15. 3 p now follows. □ 

If Q is smooth, then A is a pseudo-differential operator of order one. In this case the LP 
boundedness of the commutator [A, Xj\ in the next theorem is well known. 

Theorem 5.2. Let C = —A and Q be a bounded C 1 domain. Then 

<C p \\f\\ lp (an) (5.5) 

for any 1 < p < oo and f G L p (dQ). Lf Q is Lipschitz, the estimate 115 . 5\) holds forp = 2. 

Proof. Let Q be a bounded C l domain. Let w and u be harmonic functions in Q with 
boundary data fxi and / respectively such that (w)*, (u)* G L p (dQ). Then A(fxj) — /A(xj) — 
XiA(f) = d{w — Xiu}/dv. Let h be a harmonic function in Q such that (h)* G LP (dQ). Since 
w — XiU = on dVt and C(w — XiU) = 2 J^-, 

f d f du 

\ ^-{w-Xiu} -h\ =2\ / — -h\ < C\\f\\ L p (au) \\h\\ LP i idn) . (5.6) 

By duality this gives (15. 5p . We remark that the inequality in (I5.6P was proved in [10] for p = 2 
in Lipschitz domains. The same argument there also gives the inequality for 1 < p < oo in 
C 1 domains. Thus we may conclude that estimate (15. 5p holds for 1 < p < oo if Q is C 1 , and 
for p = 2 if Q is Lipschitz. □ 
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